Gene coexpression networks reveal a broad role for lncRNAs in inflammatory bowel disease.

The role of long noncoding RNAs (lncRNAs) in disease is incompletely understood, but their regulation of inflammation is increasingly appreciated. We addressed the extent of lncRNA involvement in inflammatory bowel disease (IBD) using biopsy-derived RNA-sequencing data from a large cohort of deeply phenotyped patients with IBD. Weighted gene correlation network analysis revealed gene modules of lncRNAs coexpressed with protein-coding genes enriched for biological pathways, correlated with epithelial and immune cell signatures, or correlated with distal colon expression. Correlation of modules with clinical features uncovered a module correlated with disease severity, with an enriched interferon response signature containing the hub lncRNA IRF1-AS1. Connecting genes to IBD-associated single nucleotide polymorphisms (SNPs) revealed an enrichment of SNP-adjacent lncRNAs in biologically relevant modules. Ulcerative colitis-specific SNPs were enriched in distal colon-related modules, suggesting that disease-specific mechanisms may result from altered lncRNA expression. The function of the IBD-associated SNP-adjacent lncRNA IRF1-AS1 was explored in human myeloid cells, and our results suggested IRF1-AS1 promoted optimal production of TNF-α, IL-6, and IL-23. A CRISPR/Cas9-mediated activation screen in THP-1 cells revealed several lncRNAs that modulated LPS-induced TNF-α responses. Overall, this study uncovered the expression patterns of lncRNAs in IBD that identify functional, disease-relevant lncRNAs.

Multiscale protein networks systematically identify aberrant protein interactions and oncogenic regulators in seven cancer types.

Global proteomic data generated by advanced mass spectrometry (MS) technologies can help bridge the gap between genome/transcriptome and functions and hold great potential in elucidating unbiased functional models of pro-tumorigenic pathways. To this end, we collected the high-throughput, whole-genome MS data and conducted integrative proteomic network analyses of 687 cases across 7 cancer types including breast carcinoma (115 tumor samples; 10,438 genes), clear cell renal carcinoma (100 tumor samples; 9,910 genes), colorectal cancer (91 tumor samples; 7,362 genes), hepatocellular carcinoma (101 tumor samples; 6,478 genes), lung adenocarcinoma (104 tumor samples; 10,967 genes), stomach adenocarcinoma (80 tumor samples; 9,268 genes), and uterine corpus endometrial carcinoma UCEC (96 tumor samples; 10,768 genes). Through the protein co-expression network analysis, we identified co-expressed protein modules enriched for differentially expressed proteins in tumor as disease-associated pathways. Comparison with the respective transcriptome network models revealed proteome-specific cancer subnetworks associated with heme metabolism, DNA repair, spliceosome, oxidative phosphorylation and several oncogenic signaling pathways. Cross-cancer comparison identified highly preserved protein modules showing robust pan-cancer interactions and identified endoplasmic reticulum-associated degradation (ERAD) and N-acetyltransferase activity as the central functional axes. We further utilized these network models to predict pan-cancer protein regulators of disease-associated pathways. The top predicted pan-cancer regulators including RSL1D1, DDX21 and SMC2, were experimentally validated in lung, colon, breast cancer and fetal kidney cells. In summary, this study has developed interpretable network models of cancer proteomes, showcasing their potential in unveiling novel oncogenic regulators, elucidating underlying mechanisms, and identifying new therapeutic targets.

DAP12 deficiency alters microglia-oligodendrocyte communication and enhances resilience against tau toxicity.

Pathogenic tau accumulation fuels neurodegeneration in Alzheimer's disease (AD). Enhancing aging brain's resilience to tau pathology would lead to novel therapeutic strategies. DAP12 (DNAX-activation protein 12) is critically involved in microglial immune responses. Previous studies have showed that mice lacking DAP12 in tauopathy mice exhibit higher tau pathology but are protected from tau-induced cognitive deficits. However, the exact mechanism remains elusive. Our current study uncovers a novel resilience mechanism via microglial interaction with oligodendrocytes. Despite higher tau inclusions, Dap12 deletion curbs tau-induced brain inflammation and ameliorates myelin and synapse loss. Specifically, removal of Dap12 abolished tau-induced disease-associated clusters in microglia (MG) and intermediate oligodendrocytes (iOli), which are spatially correlated with tau pathology in AD brains. Our study highlights the critical role of interactions between microglia and oligodendrocytes in tau toxicity and DAP12 signaling as a promising target for enhancing resilience in AD.

DAP12 deficiency alters microglia-oligodendrocyte communication and enhances resilience against tau toxicity.

Pathogenic tau accumulation fuels neurodegeneration in Alzheimer's disease (AD). Enhancing aging brain's resilience to tau pathology would lead to novel therapeutic strategies. DAP12 (DNAX-activation protein 12) is critically involved in microglial immune responses. Previous studies have showed that mice lacking DAP12 in tauopathy mice exhibit higher tau pathology but are protected from tau-induced cognitive deficits. However, the exact mechanism remains elusive. Our current study uncovers a novel resilience mechanism via microglial interaction with oligodendrocytes. Despite higher tau inclusions, Dap12 deletion curbs tau-induced brain inflammation and ameliorates myelin and synapse loss. Specifically, removal of Dap12 abolished tau-induced disease-associated clusters in microglia (MG) and intermediate oligodendrocytes (iOli), which are spatially correlated with tau pathology in AD brains. Our study highlights the critical role of interactions between microglia and oligodendrocytes in tau toxicity and DAP12 signaling as a promising target for enhancing resilience in AD.

Single-cell transcriptomic and neuropathologic analysis reveals dysregulation of the integrated stress response in progressive supranuclear palsy.

Progressive supranuclear palsy (PSP) is a sporadic neurodegenerative tauopathy variably affecting brainstem and cortical structures and characterized by tau inclusions in neurons and glia. The precise mechanism whereby these protein aggregates lead to cell death remains unclear. To investigate the contribution of these different cellular abnormalities to PSP pathogenesis, we performed single-nucleus RNA sequencing and analyzed 45,559 high quality nuclei targeting the subthalamic nucleus and adjacent structures from human post-mortem PSP brains with varying degrees of pathology compared to controls. Cell-type specific differential expression and pathway analysis identified both common and discrete changes in numerous pathways previously implicated in PSP and other neurodegenerative disorders. This included EIF2 signaling, an adaptive pathway activated in response to diverse stressors, which was the top activated pathway in vulnerable cell types. Using immunohistochemistry, we found that activated eIF2α was positively correlated with tau pathology burden in vulnerable brain regions. Multiplex immunofluorescence localized activated eIF2α positivity to hyperphosphorylated tau (p-tau) positive neurons and ALDH1L1-positive astrocytes, supporting the increased transcriptomic EIF2 activation observed in these vulnerable cell types. In conclusion, these data provide insights into cell-type-specific pathological changes in PSP and support the hypothesis that failure of adaptive stress pathways play a mechanistic role in the pathogenesis and progression of PSP.

Master Transcription Factor Reprogramming Unleashes Selective Translation Promoting Castration Resistance and Immune Evasion in Lethal Prostate Cancer.

Signaling rewiring allows tumors to survive therapy. Here we show that the decrease of the master regulator microphthalmia transcription factor (MITF) in lethal prostate cancer unleashes eukaryotic initiation factor 3B (eIF3B)-dependent translation reprogramming of key mRNAs conferring resistance to androgen deprivation therapy (ADT) and promoting immune evasion. Mechanistically, MITF represses through direct promoter binding eIF3B, which in turn regulates the translation of specific mRNAs. Genome-wide eIF3B enhanced cross-linking immunoprecipitation sequencing (eCLIP-seq) showed specialized binding to a UC-rich motif present in subsets of 5' untranslated regions. Indeed, translation of the androgen receptor and major histocompatibility complex I (MHC-I) through this motif is sensitive to eIF3B amount. Notably, pharmacologic targeting of eIF3B-dependent translation in preclinical models sensitizes prostate cancer to ADT and anti-PD-1 therapy. These findings uncover a hidden connection between transcriptional and translational rewiring promoting therapy-refractory lethal prostate cancer and provide a druggable mechanism that may transcend into effective combined therapeutic strategies. SIGNIFICANCE: Our study shows that specialized eIF3B-dependent translation of specific mRNAs released upon downregulation of the master transcription factor MITF confers castration resistance and immune evasion in lethal prostate cancer. Pharmacologic targeting of this mechanism delays castration resistance and increases immune-checkpoint efficacy. This article is featured in Selected Articles from This Issue, p. 2489.

Dysfunction of ubiquitin protein ligase MYCBP2 leads to cell resilience in human breast cancers.

Breast cancer is the most common type of cancer among women worldwide, and it is estimated that 294 000 new diagnoses and 37 000 deaths will occur each year in the United States alone by 2030. Large-scale genomic studies have identified a number of genetic loci with alterations in breast cancer. However, identification of the genes that are critical for tumorgenicity still remains a challenge. Here, we perform a comprehensive functional multi-omics analysis of somatic mutations in breast cancer and identify previously unknown key regulators of breast cancer tumorgenicity. We identify dysregulation of MYCBP2, an E3 ubiquitin ligase and an upstream regulator of mTOR signaling, is accompanied with decreased disease-free survival. We validate MYCBP2 as a key target through depletion siRNA using in vitro apoptosis assays in MCF10A, MCF7 and T47D cells. We demonstrate that MYCBP2 loss is associated with resistance to apoptosis from cisplatin-induced DNA damage and cell cycle changes, and that CHEK1 inhibition can modulate MYCBP2 activity and caspase cleavage. Furthermore, we show that MYCBP2 knockdown is associated with transcriptomic responses in TSC2 and in apoptosis genes and interleukins. Therefore, we show that MYCBP2 is an important genetic target that represents a key node regulating multiple molecular pathways in breast cancer corresponding with apparent drug resistance in our study.

Plasma proteomics of SARS-CoV-2 infection and severity reveals impact on Alzheimer's and coronary disease pathways.

Identification of proteins dysregulated by COVID-19 infection is critically important for better understanding of its pathophysiology, building prognostic models, and identifying new targets. Plasma proteomic profiling of 4,301 proteins was performed in two independent datasets and tested for the association for three COVID-19 outcomes (infection, ventilation, and death). We identified 1,449 proteins consistently associated in both datasets with any of these three outcomes. We subsequently created highly accurate models that distinctively predict infection, ventilation, and death. These proteins were enriched in specific biological processes including cytokine signaling, Alzheimer's disease, and coronary artery disease. Mendelian randomization and gene network analyses identified eight causal proteins and 141 highly connected hub proteins including 35 with known drug targets. Our findings provide distinctive prognostic biomarkers for two severe COVID-19 outcomes, reveal their relationship to Alzheimer's disease and coronary artery disease, and identify potential therapeutic targets for COVID-19 outcomes.

Multiscale network analysis identifies potential receptors for SARS-CoV-2 and reveals their tissue-specific and age-dependent expression.

The coronavirus disease 2019 (COVID-19) pandemic has affected tens of millions of individuals and caused hundreds of thousands of deaths worldwide. Here, we present a comprehensive, multiscale network analysis of the transcriptional response to the virus. In particular, we focused on key regulators, cell receptors, and host processes that were hijacked by the virus for its advantage. ACE2-controlled processes involved CD300e (a TYROBP receptor) as a key regulator and the activation of IL-2 pro-inflammatory cytokine signaling. We further investigated the age dependency of such receptors in different tissues. In summary, this study provides novel insights into the gene regulatory organization during the SARS-CoV-2 infection and the tissue-specific, age-dependent expression of the cell receptors involved in COVID-19.

Pseudo-temporal dynamics of chemoresistant triple negative breast cancer cells reveal EGFR/HER2 inhibition as synthetic lethal during mid-neoadjuvant chemotherapy.

In the absence of targetable hormonal axes, chemoresistance for triple-negative breast cancer (TNBC) often compromises patient outcomes. To investigate the underlying tumor dynamics, we performed trajectory analysis on the single-nuclei RNA-seq (snRNA-seq) of chemoresistant tumor clones during neoadjuvant chemotherapy (NAC). It revealed a common tumor trajectory across multiple patients with HER2-like expansions during NAC. Genome-wide CRISPR-Cas9 knock-out on mammary epithelial cells revealed chemosensitivity-promoting knock-outs were up-regulated along the tumor trajectory. Furthermore, we derived a consensus gene signature of TNBC chemoresistance by comparing the trajectory transcriptome with chemoresistant transcriptomes from TNBC cell lines and poor prognosis patient samples to predict FDA-approved drugs, including afatinib (pan-HER inhibitor), targeting the consensus signature. We validated the synergistic efficacy of afatinib and paclitaxel in chemoresistant TNBC cells and confirmed pharmacological suppression of the consensus signature. The study provides a dynamic model of chemoresistant tumor transcriptome, and computational framework for pharmacological intervention.

Genome-wide methylomic regulation of multiscale gene networks in Alzheimer's disease.

INTRODUCTION: Recent studies revealed the association of abnormal methylomic changes with Alzheimer's disease (AD) but there is a lack of systematic study of the impact of methylomic alterations over the molecular networks underlying AD. METHODS: We profiled genome-wide methylomic variations in the parahippocampal gyrus from 201 post mortem control, mild cognitive impaired, and AD brains. RESULTS: We identified 270 distinct differentially methylated regions (DMRs) associated with AD. We quantified the impact of these DMRs on each gene and each protein as well as gene and protein co-expression networks. DNA methylation had a profound impact on both AD-associated gene/protein modules and their key regulators. We further integrated the matched multi-omics data to show the impact of DNA methylation on chromatin accessibility, which further modulates gene and protein expression. DISCUSSION: The quantified impact of DNA methylation on gene and protein networks underlying AD identified potential upstream epigenetic regulators of AD. HIGHLIGHTS: A cohort of DNA methylation data in the parahippocampal gyrus was developed from 201 post mortem control, mild cognitive impaired, and Alzheimer's disease (AD) brains. Two hundred seventy distinct differentially methylated regions (DMRs) were found to be associated with AD compared to normal control. A metric was developed to quantify methylation impact on each gene and each protein. DNA methylation was found to have a profound impact on not only the AD-associated gene modules but also key regulators of the gene and protein networks. Key findings were validated in an independent multi-omics cohort in AD. The impact of DNA methylation on chromatin accessibility was also investigated by integrating the matched methylomic, epigenomic, transcriptomic, and proteomic data.

Multi-omics-based analysis of high grade serous ovarian cancer subtypes reveals distinct molecular processes linked to patient prognosis.

Despite advancements in treatment, high-grade serous ovarian cancer (HGSOC) is still characterized by poor patient outcomes. To understand the molecular heterogeneity of this disease, which underlies the challenge in selecting optimal treatments for HGSOC patients, we have integrated genomic, transcriptomic, and epigenetic information to identify seven new HGSOC subtypes using a multiscale clustering method. These subtypes not only have significantly distinct overall survival, but also exhibit unique patterns of gene expression, microRNA expression, DNA methylation, and copy number alterations. As determined by our analysis, patients with similar clinical outcomes have distinct profiles of activated or repressed cellular processes, including cell cycle, epithelial-to-mesenchymal transition, immune activation, interferon response, and cilium organization. Furthermore, we performed a multiscale gene co-expression network analysis to identify subtype-specific key regulators and predicted optimal targeted therapies based on subtype-specific gene expression. In summary, this study provides new insights into the cellular heterogeneity of the HGSOC genomic, epigenetic, and transcriptomic landscapes and provides a basis for future studies into precision medicine for HGSOC patients.

Biopsy and blood-based molecular biomarker of inflammation in IBD.

OBJECTIVE: IBD therapies and treatments are evolving to deeper levels of remission. Molecular measures of disease may augment current endpoints including the potential for less invasive assessments. DESIGN: Transcriptome analysis on 712 endoscopically defined inflamed (Inf) and 1778 non-inflamed (Non-Inf) intestinal biopsies (n=498 Crohn's disease, n=421 UC and 243 controls) in the Mount Sinai Crohn's and Colitis Registry were used to identify genes differentially expressed between Inf and Non-Inf biopsies and to generate a molecular inflammation score (bMIS) via gene set variance analysis. A circulating MIS (cirMIS) score, reflecting intestinal molecular inflammation, was generated using blood transcriptome data. bMIS/cirMIS was validated as indicators of intestinal inflammation in four independent IBD cohorts. RESULTS: bMIS/cirMIS was strongly associated with clinical, endoscopic and histological disease activity indices. Patients with the same histologic score of inflammation had variable bMIS scores, indicating that bMIS describes a deeper range of inflammation. In available clinical trial data sets, both scores were responsive to IBD treatment. Despite similar baseline endoscopic and histologic activity, UC patients with lower baseline bMIS levels were more likely treatment responders compared with those with higher levels. Finally, among patients with UC in endoscopic and histologic remission, those with lower bMIS levels were less likely to have a disease flare over time. CONCLUSION: Transcriptionally based scores provide an alternative objective and deeper quantification of intestinal inflammation, which could augment current clinical assessments used for disease monitoring and have potential for predicting therapeutic response and patients at higher risk of disease flares.

Plasma proteomics of SARS-CoV-2 infection and severity reveals impact on Alzheimer and coronary disease pathways.

Identification of the plasma proteomic changes of Coronavirus disease 2019 (COVID-19) is essential to understanding the pathophysiology of the disease and developing predictive models and novel therapeutics. We performed plasma deep proteomic profiling from 332 COVID-19 patients and 150 controls and pursued replication in an independent cohort (297 cases and 76 controls) to find potential biomarkers and causal proteins for three COVID-19 outcomes (infection, ventilation, and death). We identified and replicated 1,449 proteins associated with any of the three outcomes (841 for infection, 833 for ventilation, and 253 for death) that can be query on a web portal ( https://covid.proteomics.wustl.edu/ ). Using those proteins and machine learning approached we created and validated specific prediction models for ventilation (AUC>0.91), death (AUC>0.95) and either outcome (AUC>0.80). These proteins were also enriched in specific biological processes, including immune and cytokine signaling (FDR ≤ 3.72×10 -14 ), Alzheimer's disease (FDR ≤ 5.46×10 -10 ) and coronary artery disease (FDR ≤ 4.64×10 -2 ). Mendelian randomization using pQTL as instrumental variants nominated BCAT2 and GOLM1 as a causal proteins for COVID-19. Causal gene network analyses identified 141 highly connected key proteins, of which 35 have known drug targets with FDA-approved compounds. Our findings provide distinctive prognostic biomarkers for two severe COVID-19 outcomes (ventilation and death), reveal their relationship to Alzheimer's disease and coronary artery disease, and identify potential therapeutic targets for COVID-19 outcomes.

Culture shock: microglial heterogeneity, activation, and disrupted single-cell microglial networks in vitro.

BACKGROUND: Microglia, the resident immune cells of the brain, play a critical role in numerous diseases, but are a minority cell type and difficult to genetically manipulate in vivo with viral vectors and other approaches. Primary cultures allow a more controlled setting to investigate these cells, but morphological and transcriptional changes upon removal from their normal brain environment raise many caveats from in vitro studies. METHODS: To investigate whether cultured microglia recapitulate in vivo microglial signatures, we used single-cell RNA sequencing (scRNAseq) to compare microglia freshly isolated from the brain to primary microglial cultures. We performed cell population discovery, differential expression analysis, and gene co-expression module analysis to compare signatures between in vitro and in vivo microglia. We constructed causal predictive network models of transcriptional regulators from the scRNAseq data and identified a set of potential key drivers of the cultured phenotype. To validate this network analysis, we knocked down two of these key drivers, C1qc and Prdx1, in primary cultured microglia and quantified changes in microglial activation markers. RESULTS: We found that, although often assumed to be a relatively homogenous population of cells in culture, in vitro microglia are a highly heterogeneous population consisting of distinct subpopulations of cells with transcriptional profiles reminiscent of macrophages and monocytes, and are marked by transcriptional programs active in neurodegeneration and other disease states. We found that microglia in vitro presented transcriptional activation of a set of "culture shock genes" not found in freshly isolated microglia, characterized by strong upregulation of disease-associated genes including Apoe, Lyz2, and Spp1, and downregulation of homeostatic microglial markers, including Cx3cr1, P2ry12, and Tmem119. Finally, we found that cultured microglia prominently alter their transcriptional machinery modulated by key drivers from the homeostatic to activated phenotype. Knockdown of one of these drivers, C1qc, resulted in downregulation of microglial activation genes Lpl, Lyz2, and Ccl4. CONCLUSIONS: Overall, our data suggest that when removed from their in vivo home environment, microglia suffer a severe case of "culture shock", drastically modulating their transcriptional regulatory network state from homeostatic to activated through upregulation of modules of culture-specific genes. Consequently, cultured microglia behave as a disparate cell type that does not recapitulate the homeostatic signatures of microglia in vivo. Finally, our predictive network model discovered potential key drivers that may convert activated microglia back to their homeostatic state, allowing for more accurate representation of in vivo states in culture. Knockdown of key driver C1qc partially attenuated microglial activation in vitro, despite C1qc being only weakly upregulated in culture. This suggests that even genes that are not strongly differentially expressed across treatments or preparations may drive downstream transcriptional changes in culture.

Guidelines for bioinformatics of single-cell sequencing data analysis in Alzheimer's disease: review, recommendation, implementation and application.

Alzheimer's disease (AD) is the most common form of dementia, characterized by progressive cognitive impairment and neurodegeneration. Extensive clinical and genomic studies have revealed biomarkers, risk factors, pathways, and targets of AD in the past decade. However, the exact molecular basis of AD development and progression remains elusive. The emerging single-cell sequencing technology can potentially provide cell-level insights into the disease. Here we systematically review the state-of-the-art bioinformatics approaches to analyze single-cell sequencing data and their applications to AD in 14 major directions, including 1) quality control and normalization, 2) dimension reduction and feature extraction, 3) cell clustering analysis, 4) cell type inference and annotation, 5) differential expression, 6) trajectory inference, 7) copy number variation analysis, 8) integration of single-cell multi-omics, 9) epigenomic analysis, 10) gene network inference, 11) prioritization of cell subpopulations, 12) integrative analysis of human and mouse sc-RNA-seq data, 13) spatial transcriptomics, and 14) comparison of single cell AD mouse model studies and single cell human AD studies. We also address challenges in using human postmortem and mouse tissues and outline future developments in single cell sequencing data analysis. Importantly, we have implemented our recommended workflow for each major analytic direction and applied them to a large single nucleus RNA-sequencing (snRNA-seq) dataset in AD. Key analytic results are reported while the scripts and the data are shared with the research community through  GitHub. In summary, this comprehensive review provides insights into various approaches to analyze single cell sequencing data and offers specific guidelines for study design and a variety of analytic directions. The review and the accompanied software tools will serve as a valuable resource for studying cellular and molecular mechanisms of AD, other diseases, or biological systems at the single cell level.

The landscape of human tissue and cell type specific expression and co-regulation of senescence genes.

BACKGROUND: Cellular senescence is a complex stress response that impacts cellular function and organismal health. Multiple developmental and environmental factors, such as intrinsic cellular cues, radiation, oxidative stress, oncogenes, and protein accumulation, activate genes and pathways that can lead to senescence. Enormous efforts have been made to identify and characterize senescence genes (SnGs) in stress and disease systems. However, the prevalence of senescent cells in healthy human tissues and the global SnG expression signature in different cell types are poorly understood. METHODS: This study performed an integrative gene network analysis of bulk and single-cell RNA-seq data in non-diseased human tissues to investigate SnG co-expression signatures and their cell-type specificity. RESULTS: Through a comprehensive transcriptomic network analysis of 50 human tissues in the Genotype-Tissue Expression Project (GTEx) cohort, we identified SnG-enriched gene modules, characterized SnG co-expression patterns, and constructed aggregated SnG networks across primary tissues of the human body. Our network approaches identified 51 SnGs highly conserved across the human tissues, including CDKN1A (p21)-centered regulators that control cell cycle progression and the senescence-associated secretory phenotype (SASP). The SnG-enriched modules showed remarkable cell-type specificity, especially in fibroblasts, endothelial cells, and immune cells. Further analyses of single-cell RNA-seq and spatial transcriptomic data independently validated the cell-type specific SnG signatures predicted by the network analysis. CONCLUSIONS: This study systematically revealed the co-regulated organizations and cell type specificity of SnGs in major human tissues, which can serve as a blueprint for future studies to map senescent cells and their cellular interactions in human tissues.

Molecular Signature Predictive of Long-Term Liver Fibrosis Progression to Inform Antifibrotic Drug Development.

BACKGROUND & AIMS: There is a major unmet need to assess the prognostic impact of antifibrotics in clinical trials because of the slow rate of liver fibrosis progression. We aimed to develop a surrogate biomarker to predict future fibrosis progression. METHODS: A fibrosis progression signature (FPS) was defined to predict fibrosis progression within 5 years in patients with hepatitis C virus and nonalcoholic fatty liver disease (NAFLD) with no to minimal fibrosis at baseline (n = 421) and was validated in an independent NAFLD cohort (n = 78). The FPS was used to assess response to 13 candidate antifibrotics in organotypic ex vivo cultures of clinical fibrotic liver tissues (n = 78) and cenicriviroc in patients with nonalcoholic steatohepatitis enrolled in a clinical trial (n = 19, NCT02217475). A serum protein-based surrogate FPS was developed and tested in a cohort of compensated cirrhosis patients (n = 122). RESULTS: A 20-gene FPS was defined and validated in an independent NAFLD cohort (adjusted odds ratio, 10.93; area under the receiver operating characteristic curve, 0.86). Among computationally inferred fibrosis-driving FPS genes, BCL2 was confirmed as a potential pharmacologic target using clinical liver tissues. Systematic ex vivo evaluation of 13 candidate antifibrotics identified rational combination therapies based on epigallocatechin gallate, which were validated for enhanced antifibrotic effect in ex vivo culture of clinical liver tissues. In patients with nonalcoholic steatohepatitis treated with cenicriviroc, FPS modulation was associated with 1-year fibrosis improvement accompanied by suppression of the E2F pathway. Induction of the PPARα pathway was absent in patients without fibrosis improvement, suggesting a benefit of combining PPARα agonism to improve the antifibrotic efficacy of cenicriviroc. A 7-protein serum protein-based surrogate FPS was associated with the development of decompensation in cirrhosis patients. CONCLUSION: The FPS predicts long-term fibrosis progression in an etiology-agnostic manner, which can inform antifibrotic drug development.

Integrative metabolomics-genomics approach reveals key metabolic pathways and regulators of Alzheimer's disease.

Metabolites, the biochemical products of the cellular process, can be used to measure alterations in biochemical pathways related to the pathogenesis of Alzheimer's disease (AD). However, the relationships between systemic abnormalities in metabolism and the pathogenesis of AD are poorly understood. In this study, we aim to identify AD-specific metabolomic changes and their potential upstream genetic and transcriptional regulators through an integrative systems biology framework for analyzing genetic, transcriptomic, metabolomic, and proteomic data in AD. Metabolite co-expression network analysis of the blood metabolomic data in the Alzheimer's Disease Neuroimaging Initiative (ADNI) shows short-chain acylcarnitines/amino acids and medium/long-chain acylcarnitines are most associated with AD clinical outcomes, including episodic memory scores and disease severity. Integration of the gene expression data in both the blood from the ADNI and the brain from the Accelerating Medicines Partnership Alzheimer's Disease (AMP-AD) program reveals ABCA1 and CPT1A are involved in the regulation of acylcarnitines and amino acids in AD. Gene co-expression network analysis of the AMP-AD brain RNA-seq data suggests the CPT1A- and ABCA1-centered subnetworks are associated with neuronal system and immune response, respectively. Increased ABCA1 gene expression and adiponectin protein, a regulator of ABCA1, correspond to decreased short-chain acylcarnitines and amines in AD in the ADNI. In summary, our integrated analysis of large-scale multiomics data in AD systematically identifies novel metabolites and their potential regulators in AD and the findings pave a way for not only developing sensitive and specific diagnostic biomarkers for AD but also identifying novel molecular mechanisms of AD pathogenesis.

DDR2 coordinates EMT and metabolic reprogramming as a shared effector of FOXQ1 and SNAI1.

While multiple transcription factors (TFs) have been recognized to drive epithelial-mesenchymal transition (EMT) in cancer, their interdependence and context-dependent functions are poorly understood. In this study, we show that FOXQ1 and SNAI1 act as independent TFs within the EMT program with a shared ability to upregulate common EMT TFs without reciprocally impacting the expression of one another. Despite this independence, human mammary epithelial cells (HMLE) with ectopic expression of either FOXQ1 or SNAI1 share a common gene set that is enriched for a DDR2 coexpression signature. Further analysis identified DDR2 as the most upregulated receptor tyrosine kinase and a shared downstream effector of FOXQ1 and SNAI1 in triple-negative breast cancer (TNBC) cell lines. Alteration of DDR2 expression in either FOXQ1 or SNAI1 driven EMT models or in TNBC cells resulted in a profound change of cell motility without significantly impacting EMT marker expression, cell morphology, or the stem cell population. Lastly, we demonstrated that knockdown of DDR2 in the FOXQ1-driven EMT model and TNBC cell line significantly altered the global metabolic profile, including glutamine-glutamate and Aspartic acid recycling.